Twirling motion of actin filaments in gliding assays with non-processive myosin motors 
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We present a model study of gliding assays in which actin filaments are moved by non-processive 
myosin motors. We show that even if the power stroke of the motor protein has no lateral component, 
the filaments will rotate around their axis while moving over the surface. Notably, the handedness 
of this twirling motion is opposite from that of the actin filament structure. It stems from the fact 
that the gliding actin filament has "target zones" where its subunits point towards the surface and 
are therefore more accessible for myosin heads. Each myosin head has a higher binding probability 
before it reaches the center of the target zone than afterwards, which results in a left-handed twirling. 
We present a stochastic simulation and an approximative analytical solution. The calculated pitch 
of the twirling motion depends on the filament velocity (ATP concentration). It reaches about 
400nm for low speeds and increases with higher speeds. 



Introduction 

Gliding assays, also known as motility assays, repre- 
sent the oldest in vitro technique to study motor proteins 
P, Q • They consist of attaching motors (like myosins or 
kinesins) with their tails to a glass surface and adding the 
filaments (actin or microtubules). The motors will then 
pull the filaments and make them glide over the surface 
(Fig. [T]A_). Gliding assays are the most convenient way 
of testing motors for their functionality, measuring their 
speed in the absence of load and for testing their pro- 
cessivity. Several experimental and theoretical studies 
were dealing with the pathways of such filaments in the 
two-dimensional plane [3,, i^, ^5,] . Interestingly, one group 
observed that gliding actin filaments move in a helical 
fashion f^, 0| • In a subsequent experiment the pitch of 
rotation was determined as about 1 /xm, although the ap- 
plied optical detection method did not allow discrimina- 
tion between left- and right-handed rotation Q . 

Helical motion of myosin motors has been very im- 
portant in a somewhat different context. The processive 
motor myosin V has an average step size that is close, 
but not precisely equal to the actin periodicity. The heli- 
cal motion of a motor around the actin filament therefore 
presents a very accurate way of measuring the difference 
between its step size and the filament pitch. Ali and 
coworkers have observed that myosin V walks on an actin 
filament along a left-handed helix with a pitch of 2.2 /im 
and thus has a step size slightly shorter than the actin 
half-pitch (for a discussion of the myosin V step size see 

Tio, 11]). 

Myosin VI, despite having a shorter lever arm than 
myosin V, showed either straight walking, or, in 20% of 
cases, a helical path with a pitch of 2.3 /im 12"]. Sun 
et al. [l^ confirmed this result, but also showed that 
the relatively straight motion contains a large amount 
of random wiggling. New experiments on myosin X also 
show a left-handed helical motion with a pitch that is 
somewhat shorter than that of myosin V and VI [l3| . 

In a recent experimental study Beausang and cowork- 
ers [iH] used polarized total internal reflection microscopy 
to study the twirling motion of actin filaments in gliding 



assays with processive myosin V and non-processive mus- 
cle myosin (myosin II). While the twirling of filaments 
driven by myosin V agreed with the helical movement 
of single molecules mentioned above, myosin II interest- 
ingly showed a left-handed twirling motion as well. This 
result came as a surprise and the left-handed rotation is 
opposite from the observations by Nishizaka et al. Q. 
But they are not in direct contradiction, as they were 
obtained with quite different ATP concentrations. 

While the pitch of the twirling motion is a direct mea- 
sure for the step size of processive motors, its interpreta- 
tion is more complicated with non-processive ones. They 
could clearly generate twirling motion if there was a lat- 
eral component of the power stroke. In fact, there exists 
indirect evidence for such an asymmetry in myosin V . 
However, we will show in this paper that there is another, 
more subtle, effect that can cause twirling motion of actin 
filaments in a gliding assay, even if the myosin heads ex- 
hibit no lateral motion. This effect stems from the fact 
that myosin heads can only bind to an actin filament in 
so called target zones, where the actin binding sites have 
approximately the right orientation (Fig. [Tp) [l3, [l^ . 
When a target zone is approaching a myosin head, the 
latter is more likely to bind at the beginning of the target 
zone than at its end, because it is more likely that it is 
already bound by that time. In this paper we will show 
simulation results and develop an approximative theory 
to estimate the pitch of helical motion resulting from this 
effect. Of course, we cannot exclude that there are other 
contributions towards the helicity. But, because the rota- 
tion is relatively weak as compared with the longitudinal 
motion, these effects can easily be treated separately and 
the total rotation is simply their superposition. 



Model definition 

In order to concentrate on the effect of filament ro- 
tation, we define a model with simplified myosin kinet- 
ics, essentially containing a detached state, a bound pre- 
powerstroke state, a bound post-powerstroke state with 
ADP and a bound rigor state. At the same time, we take 
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FIG. 1: (A) Model definition. The linear motion and rota- 
tion of the actin filament are denoted with X and O. Motors 
are situated underneath the actin filament at randomly dis- 
tributed positions x-'j^j and each motor can be in the detached, 
pre-powerstroke or post-powerstroke state. (B) Duty cycle of 
a myosin head. The head binds to the actin site i with rate k\, 
then quickly undergoes the power stroke with rate kps, re- 
leases ADP with rate k-ADP and detaches with rate ko- The 
step size is d = 8 nm. (C) Front view of a bound head. For the 
sake of simplicity, we assume that the azimuthal component 
of the elastic energy needed to bind to site i only depends on 
the angle i9. (D) Target zones on an actin filament. When a 
filament is moving past an actin head, it is more likely to bind 
to submits pointing towards the myosin-covered surface. Red 
color denotes the sites with the highest binding rate. Besides 
that, sites at the beginning of a target zone have a higher 
binding probability than those at its end. 



into account the full helical structure of the actin fila- 
ment. We propose an actin filament which can move in 
one direction and rotate around its axis. The position of 
the actin filament at a given time is therefore described 
with the coordinates {X, 8). As follows from the helical 
structure of the actin filament, the x-coordinate of each 
bindings site is X + ia and its azimuth angle Q + ido, 
where do = -(13/28) x 360° = -167.14° is the rotation 
and a = 2.75 nm the axial rise per subunit. The defini- 
tion of the model is illustrated in Fig. [1] 

We assume that the myosin motors are distributed ran- 
domly directly under the gliding actin filament (a discus- 
sion how this simplified, one-dimensional model follows 
from the full, 2-D model is given in the Appendix). The 
motor numbered j is anchored at position x'^j^. The elas- 
tic energy cost of binding a head to the site i consists 
of a longitudinal component with stiffness K and and an 
angular component with stiffness and can be written 
as 



1 



K{X 



^Ai? + lK4e + ii9o + 27rnf (1) 



with n chosen such that the angle 8 + i-do + 2-Kn falls 
into the interval [— tt, tt]. The binding rate is then pro- 
portional to the Boltzmann factor 



kA exp 



(2) 



This is essentially the expression used by Steffen et al. 
[l3| to fit binding rates of a single myosin head to the 
actin filament. They determined the value of the angular 
stiffness expressed with the dimensionless coefficient 



a — K^/ksT 



(3) 



as a = 3.7. However, this value needs to be regarded 
as a lower estimate, as it might partially result from 
torsional compliance of the actin filament, rather than 
myosin heads. We therefore use three different values of 
a = 4, 6 and 8 in the simulation. The angular contribu- 
tion to the Boltzmann factor for a set of binding sites is 
shown in Fig. [2l For the longitudinal compliance, we use 
the value K ~ 0.5pN/nm, somewhat below the stiffness 
of myosin heads in muscle, which is about 2.5pN/nm 
p^ . The lower stiffness reflects the additional compli- 
ance due to myosin tails and roughly corresponds to the 
value obtained with optical tweezers, 0.69 ± 0.47pN/nm 
0. 

The force and the torque that a myosin head numbered 
j, bound to site i, exerts on the filament are 

= K{x{j + S-X-ia) = -Ki}{Q + ido + 2Tm) . 

(4) 

Here we introduced the displacement 5 that has the value 
in the pre-powerstroke and d = 8 nm in the post- 
powerstroke and rigor state. 

The simplified model for the duty cycle of the myosin 
head is defined as follows. A head binds to an actin site 
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FIG. 2: Target zones on an actin filament. The diagram 
shows the binding rate to the particular site if there is no lon- 
gitudinal strain, i.e., if the head is horizontally aligned with 
the binding site. The binding rate has its maxima for actin 
subunits oriented downwards (towards the myosin-coated sur- 
face). The dimensionless angular stiffness is q = 4. 



with the rate ki given by Eq. The power stroke (a 
transition from S — do S — d) takes place with the 
rate kpg . It is followed by the release of ADP with the 
rate k-ADP- Detachment follows after binding a new 
ATP molecule, therefore its rate depends on the ATP 
concentration, kn — fc£)[ATP]. We neglect the strain de- 
pendence of those rates, as well as the existence of the 
reverse transitions. It should be noted that in this for- 
mulation the model is not thermodynamically consistent. 
However, as we are only interested in dynamics at low 
loads, this does not significantly affect the results. 

We assume that the filament position is quickly equi- 
librated after each step, therefore it always fulfills 
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(5) 



Simulation results 



Given the known structure of the actin helix and the 
power stroke size of myosin, our model essentially has two 
important parameters: the angular stiffness of myosin 
heads, K^, and the ratio between the detachment- and 
the attachment rate, ko/kA- The latter is a function of 
the ATP concentration and is closely related to the duty 
ratio of motors. For other parameters we use the val- 
ues given in Table [H The power stroke, connected with 
the phosphate (Pi) release, is assigned a very fast rate, 
kps = 10000 s~^, and can be considered as taking place 
immediately after binding. The maximum attachment 
rate kA, i.e., the attachment rate to sites that do not re- 
quire any elastic distortion, can be estimated as 50 s~^. 
This reflects the estimated average attachment rate of 
30 s~^, or a maximum ATP turnover rate of 25 s~^ in 
muscle [2l|. The ADP release takes place with the rate 
k-ADP = 1000 s^^, characteristic for the fast myosin iso- 
form For the detachment rate, which is determined 
by the ATP binding rate, we use fc^ = 5 fiM'^s'^ p^ . 



kA 50 s~^ Attachment rate 

kps 10000 s~^ Power stroke rate 

k-ADP 1000 s~'^ ADP release rate 

ko (5 AtM~^s~^) [ATP] Detachment rate 

d 8nm Power stroke size 

K 0.5pN/nm Myosin stiffness (longitudinal) 

Ki) — akpT Myosin stiffness (angular) 

a 4, 6, 8 Dimensionless angular stiffness 

a 2.75 nm Distance between actin subunits 

i9o —167.14° Angle between actin subunits 

I 5.5 nm Actin filament length 

p 20 fim~^ 1-d myosin density on surface 



fcflT 4.14 X 10" 



Thermal energy 



TABLE I: Simulation parameters 



The stochastic simulation essentially followed the fol- 
lowing algorithm: 

1. Distribute the positions of myosin motors ran- 
domly along the distance covered by the actin fila- 
ment, with an average linear density p. Set X — 
and 9 = 0. 

2. Determine the total rate of all possible transitions 
as 



k, 



total 



E 

jGMotors 



EjGBinding sites ^a(j) if ^0*01' J dctachcd 

kps if motor j in pre-PS state 

k-ADP if motor j in post-PS state 

ko if motor j in rigor state 

k\{j) is determined using Eq. ^ with the current 
values of X and 0. 

3. Determine the time until the next step as Ai — 
'^totai 1^^(1/^)1 where r is a random number between 
and 1. 

4. Choose randomly one of the possible steps (attach- 
ment, detachment, power stroke, ADP release), so 
that the probability of choosing a certain step is 
given by its rate, divided by fctotai- 

5. Change the state of the chosen motor and update 
the filament position X and angle according to 

Eq. dni). 

6. Continue with step 2 until X > Amax (Amax = 
1000 /xm). 

7. Determine the average speed as u = X/t and pitch 
as A = 27rA/e. 

The results of this numerical simulation are shown in 
Fig. [21 which shows the inverse pitch of twirling as a func- 
tion of the ATP concentration for three different values of 
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FIG. 3: (A) Inverse twirling pitch (number of rotations per 
micron traveled) as a function of the ATP concentration. The 
lines represent 3 different values of the angular stiffness {a = 
4, a = 6 and a = 8). Negative signs denote left-handed 
rotation. The upper scale shows the velocity (for a — 4). (B) 
Gliding velocity as a function of the detachment rate. (C) 
Duty ratio rj, giving the average fraction of myosin heads in 
the bound state. 




FIG. 4: Simplified model used for the analytical solution. The 
filament travels with a stationary velocity v and rotates with 
an angular velocity lj. x denotes the position of a myosin 
anchoring point relative to the center of the target zone. 



the angular stiffness K^. For reference, velocity (Fig.[5}3) 
and duty ratio (Fig.[3p) are included as well. The behav- 
ior of the pitch is non-monotonous: it has a minimum of 
about 400 — 500 nm at intermediate speeds, but increases 
both at high, as well as very low speeds. 

These results show that the helicity of the actin fil- 
ament is sufficient to explain the twirling motion in a 
gliding assay. Somewhat counter-intuitively, this rotation 
is left-handed, and therefore opposite from the handed- 
ness of the actin filament. The effect becomes weaker for 
high speeds (where the distance traveled between two at- 
tachment events becomes longer), as well as for very low 
speeds, where motors have enough time to bind even to 
unfavorable sites outside the target zones. 



Analytical approximation 

In the following we will describe an approximative ana- 
lytical solution with the aim of understanding and quan- 
titatively reproducing the twirling dynamics. The essen- 
tial simplification we will make is to neglect the discrete 
nature of binding sites on the actin filament and replace 
them with a continuous helical "groove" . The simplified 
model is shown in Fig. |4l We denote each head with its 
root position x relative to the center of the target zone: 



xm 



X 



-e + nL 



(6) 



with n such that —^<x<^.A bound head is addition- 
ally characterized by the strain ^, which is the position of 
its root relative to its binding site. For for a head bound 
to site i, this is £,i = xm — X — ia. 

In the original model the total binding rate for a motor 
positioned at x is 



kA{x) ^^k\{x) kA exp 



2kBT 



(7) 



In the sum, we only consider sites that are turned to- 
wards the motor-covered surface, therefore we can write 
the angle as 



(8) 
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L period (half-pitch) of the actin superhelix 

X motor root position relative to the center of the 

nearest target zone 
^ motor root position relative to its binding site on 

actin 

k\ attachment rate to site i 

ItA (x) total attachment rate for a motor positioned at x 
{£,a{x)) average strain of newly attached motors posi- 
tioned at X 

{£,a) average strain of all newly attached motors 
(xa) average position of newly attached motors, rela- 
tive to the target zone 
(iSa) average angular strain at attachment 
V filament velocity 

c apparent velocity of the actin helix 

u! angular velocity of actin rotation 

TABLE II: Variables used in the analytical calculation 



and by completing the square in the numerator we obtain 



We can now set up a Master equation for the prob- 
ability that a motor positioned at x is in the attached 
state 

dtA{x, t) - cd^A{x, t) ^ kAix){l ~ A{x, t)) - koAix, t) 

(13) 

and set dtA = to obtain the stationary solution. A also 
has to fulfill the periodic boundary condition 



A{-L/2) = A{LI2) 



(14) 



The expectation value of the attachment position can be 
calculated as 



/-L/2 ^kA{x){\ — A{x))dx 



J%kA{x){l - A{x))dx 



(15) 



and the average strain at the time of attachment follows 
from pTjl : 



{U) = 



K + K: 



(16) 



k\{x) = kA{x,i,i) 
kA exp 



K 



-K' -^l K+K', 



2kBT 



(9) 



In this equation we introduced the reduced angular stiff- 
ness K'^ = {-k"^ / L'^)K^. When we neglect the discreteness 
of binding sites and extend the summation beyond one 
period, we can replace the sum by an integral 



kA{x) 




kAix,^)d^ 



2'KkBT 



K + K: 



7 exp 



KK'^/{K + K'^)x^ 



2kBT 



(10) 



For a fixed x, the the expected value of the strain at the 
time of attachment can be calculated using ([9]) : 



T.^^a{^) J^^kA{x,Od^ K + K', 



•i9 

11) 



The expected value of the azimuthal angle at time of 
attachment follows from Eq. ([5]) 



K 



K + K 



r(i) 



(12) 



In the stationary state, the filament moves with ve- 
locity X = V and rotates with angular velocity Q = uj. 
From Eq. ^ it follows that x = — (v — ^uj) — — c. This 
is the apparent velocity with which the helix moves along 
the surface. 



Because the strain on a motor changes with time as ^ = 
—V, the average strain of all bound motors is {S,a) —v /ko- 
The force per motor is then F = K + d) = K{{^a} — 
v/ko + d. As the total force produced by all motors has 
to be zero, we obtain an expression for the velocity 



V = kD{d+{U)) ■ 



(17) 



The same type of calculation as for the velocity can be 
made for the angular velocity. Motors attach with an 
average angle {"Oa)- As the filament rotates, their an- 
gle changes &s ■& — lu. The average angle of all motors 
is (i^a) + iojkY) and needs to be zero because of torque 
balance, therefore 



= -k 



D 



(^.>--^.^0(..). (18) 



These equations, together with the Master equation (|13p . 
the periodic boundary condition (|14p . the expression for 
{£,a) (fT6| and for ItA Q allow us to numerically deter- 
mine the velocity v and the distribution of attached heads 
A{x) in a self-consistent manner. An example of the so- 
lution A{x), along with the attachment rate kA{x) and 
attachment flux kA{x){l—A{x)) is shown in Fig.[Sl\. Well 
visible is the asymmetry in the attachment flux. The ex- 
pectation value of x at the time of attachment ({xa}) as a 
function of the ratio ko/kA'^s shown in Fig. [S]B . It reaches 
its maximum when ko/kA is such that each motor trav- 
els an average path of ~ L/3 between two attachment 
events. 

This finally gives us the expression for the twirling 
pitch 



2'KV 



OJ 



2Ld I 



K 



{XA 



2L§ . (19) 
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FIG. 5; (A) Top panel: probability A{x) that a motor with a 
position X relative to the center of a target zone is in the bound 
state. The parameters are ko/kA = 0.01,0.1,1,10 (each a 
different hne type) and a = 4. Middle panel: dimensionless 
attachment rate kA{x)/kA- Bottom panel: attachment flux, 
obtained as the product of the probability that a motor is 
in the detached state and the attachment rate. Note that 
the asymmetry is most pronounced for ko/kA = 0.1. (B) 
The expectation value of x at the time of binding, (xa), as a 
function of ko/kA- 

The results are shown in Fig. [H] and compared with 
simulation data from Fig. [31 The simulation results are 
well reproduced, although there is a certain discrepancy 
which is more pronounced for low values of the angular 
stiffness a. The main reason for this discrepancy is the 
extrapolation beyond the boundaries of one period, which 
was used in the derivation of Eq. PU)) . Other (minor) 
sources of deviation are the neglected discrete nature of 
binding sites and of the fact that each binding site can 
only be occupied by one head at a time. 

Discussion 

In this study we have demonstrated that the helical 
actin structure, along with the fact that myosin heads 



FIG. 6: Inverse pitch — A~^, obtained from Eq. (|19p (con- 
tinuous line), and from simulation described in the previous 
section, however assuming k-ADP = oo. (circles). The top 
panel shows data for the angular stiffness a — 4, the middle 
a = 6 and the bottom a = 8. The minor deviation is mainly 
due to the extrapolations we made in regions between tar- 
get zones. The agreement is better for large a values, where 
target zones become more localized. 



preferentially bind to those sites oriented towards them, 
is sufficient to explain left-handed rotation in a gliding 
assay. The maximum twirling motion is achieved at rela- 
tively low speeds (below lOOnm/s). Twirhng is reduced 
with higher speeds, achieved at higher ATP concentra- 
tions. Interestingly, it is also reduced under extremely 
low ATP concentrations, when the velocity drops un- 
der lOnm/s. However, in this regime the results de- 
pend strongly on the choice of the angular stiffness a, 
which is less well known. The minimum pitch resulting 
from this effect lies in the order of 400-500 nm, which 
is in good agreement with recent experimental results 
[TH - The model also makes a testable prediction that 
the pitch should increase with a higher ATP concentra- 
tion. Because pitch only depends on the ratio between 
the attachment and detachment rate, addition of ADP 
should have the same effect on the pitch as a reduced 
ATP concentration that yields the same filament speed. 
If this turns out not to be the case, it will be a strong in- 
dication of a lateral conformational change in the myosin 
head connected with the release of ADP. 

Although the qualitative aspects of our theory are 
generic and practically independent of any assumptions 
other than the helical actin structure, there are alterna- 
tive effects that could well contribute to the twirling mo- 
tion. One such possibility is that the power-stroke of the 
myosin head contains a lateral (azimuthal, off-axis) com- 
ponent. A similar effect could also result from an asym- 
metric attachment rate, which could cause the attached 
motors to exert a certain torque on the filament imme- 
diately after binding. Such a torque does not contradict 
the laws of thermodynamics because the first bound state 
is not in equilibrium with the detached state. A related 
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idea is described by Beausang et al. [l5| as the "rigor drag 
model" in which heads in the rigor state exert a torque in 
the opposite direction from that immediately after bind- 
ing. This results in a pitch that depends on the fraction 
of time spent in the rigor state. The important difference 
between the two concepts is that in our model the torque 
generated by newly attached myosin heads depends on 
the ATP concentration, whereas in the rigor drag model 
this torque is constant and the variable pitch is caused 
by different dwell times in different states. 

In any case both the effect described here and the ex- 
plicit lateral component of the power stroke will be su- 
perimposed. So it is theoretically possible, even though 
the proposition is purely speculative at the moment, that 
the power stroke might have the opposite helicity, i.e., it 
would lead to a right-handed filament rotation. In such 
a case, there could be a cross-over form right handed 
motion under high ATP concentrations, to left-handed 
under low. This is could be one possibility to reconcile 
the recent results (l^ with those by Nishizaka et al. Q- 

Recent experiments also revealed rotation of micro- 
tubules moved by monomeric kinesin-1 and Eg5 [23| . 
The theory we presented in this paper is not applicable to 
microtubules, because they have no distinct target zones. 
Any rotation resulting from an effect of the kind we de- 
scribe here would be negligible. Therefore, as suggested 
in (2^ , the rotation caused by kinesins has to result from 
a lateral (off-axis) component of a power stroke, or from 
an asymmetry in the binding rate. 



Appendix 

In the following we will discuss how the simplified 
model, which assumes that motors are distributed one- 
dimensionally underneath the actin filament, and which 
we use throughout the main text, relates to the full two- 
dimensional model. In the 2-D model, which describes 
the actual situation in a gliding assay, myosin heads are 
distributed all over the glass surface and their positions 
are described with two coordinates, {x\.j , y''j^,j) . The po- 
sition of the actin filament is described with the coor- 
dinates {X, Y) and the angle G. We assume that the 
filament keeps its direction parallel to the AT-axis. The 
elastic distortion of the head j binding to the site i can 
then be written as 



The total binding rate to site i of all motors located at 
longitudinal position xm is then 



1 



C/, = ^K{X + ia- x\,Y + Ua{Q + i^a,Y - y\j) (20) 

where Ua is an unknown function of the azimuth angle of 
binding site i and of the lateral position of the filament 
relative to the motor. 



kAixu) 

kA exp 



Ui{X + ia- XM, O -I- ido, Y - i/m) 



X pD{xM,VM)dyM (21) 

where p is the 2-D surface density of myosin motors 
and D{xm tVm) the probability that a motor at that 
position is in the detached state. If we assume that 
D{xM tUm) — D{xM,'2Y — j/m), i.e., that the distribu- 
tion of unbound heads is symmetric with respect to the 
filament, the resulting function kA^XM) has to be sym- 
metric in -|- i'do. We can therefore approximate it with 
the expression used in Eqs. (|ll2p . 

The average torque generated by a head that binds to 
site i can be calculated as 



(Af) = - 



f^^ exp i-UjkBT) D{xM,yM) jdUjde) dyM 

jToo <5XP {~Ui/kBT) D{xM,yM)dyM 
ksT dkAixAi) 



kA{xM) dQ 



(22) 



This second expression is equivalent to that in the 1-D 
model. 

The average force generated by a head that binds to 
site i is determined the same way: 



(Fl) 



J^^ exp i-UjkBT) DixM,yM) jdUjdY) dyM 

jToo ^^P {-Ui/kBT) D{xM,yAi)dyM 

(23) 

If D(xM,yM) is independent of yM, or, more generally, 
if it has a dependence that can be written as a function 
of C/i, the integral in the numerator is and there is 
no lateral force. However, with different distributions 
D{xM,yM), a small lateral force is possible, so that the 
filament could show some sideways motion in the 2-D 
model. While the torque results from a D{xm, yu) which 
is asymmetric in xm, a lateral force needs asymmetry in 
both coordinates and is therefore a higher-order effect. 
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